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Abstract 

Conversion of slope cropland to perennial vegetation has a significant impact on soil organic carbon (SOC) stock in A 
horizon. However, the impact on SOC and its fraction stratification is still poorly understood in deep soil in Loess Hilly 
Region (LHR) of China. Samples were collected from three typical conversion lands, Roblnia psendoacada (RP), Caragana 
Korshinskii Kom (CK), and abandoned land (AB), which have been converted from slope croplands (SC) for 30 years in LHR. 
Contents of SOC, total nitrogen (TN), particulate organic carbon (POC), and labile organic carbon (LOC), and their 
stratification ratios (SR) and carbon management indexes (CMI) were determined on soil profiles from 0 to 200 cm. Results 
showed that the SOC, TN, POC and LOC stocks of RP were significantly higher than that of SC in soil layers of 0-10, 10-40, 
40-100 and 100-200 cm (P<0.05). Soil layer of 100-200 cm accounted for 27.38-36.62%, 25.10-32.91%, 21.59-31.69% and 
21 .08-26.83% to SOC, TN, POC and LOC stocks in lands of RP, CK and AB. SR values were >2.0 in most cases of RP, CK and 
AB. Moreover, CMI values of RP, CK, and AB increased by 11.61-61.53% in soil layer of 100-200 cm compared with SC. 
Significant positive correlations between SOC stocks and CMI or SR values of both surface soil and deep soil layers indicated 
that they were suitable indicators for soil quality and carbon changes evaluation. The Grain-to-Green Program (GTGP) had 
strong influence on improving quantity and activity of SOC pool through all soil layers of converted lands, and deep soil 
organic carbon should be considered in C cycle induced by GTGP. It was concluded that converting slope croplands to RP 
forestlands was the most efficient way for sequestering C in LHR soils. 
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Introduction 

Soil organic carbon (SOC) is a dynamic component of the 
terrestrial system, with internal changes in both vertical and 
horizontal directions and external exchanges between the atmo- 
sphere and the biosphere [1]. SOC storage is estimated at 
approximately 1.500 Pg globally, which is about two and three 
times the size of carbon pools in the atmosphere and vegetation, 
respectively [2] . Since carbon uptake and storage is tightly linked 
to the nitrogen (N) cycle, it is equally important to understand how 
N pools and fluxes are affected by land use change [3] . Moreover, 
more than 50% of the total SOC is stored in the subsoil [4]. The 
proportion of soil organic matter (SOM) stored in the first meter of 
the world soils below 30 cm depth ranges 46%~63%, except for 
Podzoluvisols, where 30% of SOC is stored below the depth of 
30 cm [4]. A recent study also suggests that in the northern 
circumpolar permafrost region, at least 6 1 % of the total soil C is 
stored below 30 cm [5] . Therefore, subsoil C may be even more 
important in terms of source or sink for COj than topsoil C [6]. 
Considering the potential role of SOC in atmospheric CO2 sink, it 
is important to understand what leads to sequestration of large 
amounts of SOC in the subsoil or even in deep soil. However, the 



SOC contents in deep soil layers are not fully understood in LHR 
of China to date. 

As an indicator of soil quality, SOM stratification, which is 
related to the rate and amount of SOC sequestration [7], is 
common in many natural ecosystems [8] and managed grasslands 
and forests [9-10]. Stratification ratio (SR) is defined as the ratio of 
a soil property at the surface layer to that at a deeper layer. In 
general, high SR values indicate good soil quality and are usually 
used to assess agricultural practices [7] . For instance, SR values for 
SOC at depths of 0-5 cm and 20-40 cm range from 1.1 to 1.5 
under traditional tillage (TT) while from 1.6 to 2.6 under 
conservation tillage (CT) [11]. Little information is available on 
natural ecosystems and managed shrubs or forests land. Addition- 
ally, under semiarid climate, SOC in active fractions, such as 
particulate organic carbon (POC) and labile organic carbon 
(LOC) was more sensitive to soil management practices than total 
SOC [12]. Previous researches have indicated that changing rate 
of POC and LOC was faster than SOC in whole soil [11], and 
they could be an early indicator for SOC change in soil [13]. 
Meanwhile, the carbon management index (CMI), which is 
derived from the total soil organic C pool and C lability, had been 
extensively used as a sensitive indicator of SOC variation rate in 
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response to soil management changes [14-15]. Therefore, under 
semi-arid climate, using SR of total SOC and of different SOC 
fractions may be useful to reveal how soil management affects soil 
quality and helpful to understand the mechanism of SOC 
transformation and cycling in subsoil as well as in deep soil. In 
LHR of China, soil erosion and desertification are causing a loss of 
net primary productivity that was estimated as high as 12 kg C 
ha 'y ' [16]. To counteract soil erosion and other environmental 
problems, an environmental protection policy was implemented 
by Chinese central government, which was known as the Grain to 
Green Program (GTGP). The purpose of GTGP was to convert 
up to 26.87 million ha low-yield sloped croplands (>25°) into 
forests, shrubs or grasslands by the end of 2008 [17]. It is the first 
and the most ambitious "payment-for-ecosystem-services" pro- 
gram in China to date [18]. Although the initial goal of GTGP was 
to control soil erosion in China, it also plays a significant role in 
circulation of SOC and total nitrogen (TN). In recent years, a few 
studies estimated the effects of GTGP on vegetation structure, 
economic benefits, soil physiochemical properties, and niche 
characteristics [19-21]. However, SR values of SOC and/or TN 
and CMI value among different land use types are rarely reported. 
Especially, information on dynamics of C in deep soil is largely 
ignored in this region. 

This study aimed to: 1) analyze the contents of SOC, TN, POC 
and LOG and their vertical distributions at the depths of 0- 
200 cm; 2) assess the stocks of SOC and TN at different soil depths 
of three land use types; and 3) evaluate the soil quality of different 
land use types using SR and CMI values as the main assessment 
parameters. 

Materials and Methods 

AU sites in the watershed we were selected for study was 
determined through interviews with local farmers (Mr. Yibin 
Zhang, Soil and Water Conservation Experiment Station, 
Northwest A&F University, Ansai County, Shaanxi, NW Chi- 
na) .We state clearly that no specific permissions were required for 
the location. We confirm that the location is not privately-owned 
or protected in any way. We confirm that the field studies do not 
involve endangered or protected species. 

Research area 

The study was conducted in the Zhifanggou catchment 
(36°46'42"-36°46'28"N, 109°13'46"-109°16'03"E), which is lo- 
cated in Ansai county, central LHR (see Fig. 1). Ansai is a typical 
county characterized by semi-arid climate and hUly loess 
landscape in the Loess Plateau with an annual average temper- 
ature of 8.8°C, and an average annual precipitation of 505 mm. 
60% of precipitation occurs between July and September 
(~300 mm in dry years while >700 mm in wet years). Accumu- 
lated temperatures above 0°C and 10°C are 3733°C and 3283°C, 
respectively. On average, there are about 157 frost-free days and 
2415 h svmny time each year. Arable farming mostly occurs on 

Table 1. Characteristics of different vegetation types. 
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Figure 1. Location of the Loess Plateau and the study site. 

dol:10.1371/joumal.pone.0099657.g001 

sloping lands without irrigation. The loess parent material at the 
site has an average thickness of approximately 50-80 m and the 
soil in this region is classified as Calciustepts soil [22]. Sand (2- 
0.05 mm) and silt (0.05-0.002 mm) account for approximately 
29.22% and 63.56% in soil depth of 0-20 cm, respectively. The 
soil is highly erodible, with an erosion modulus of 10,000— 
12,000 Mg-km ^-yr ' before the start of restoration efforts [23]. 
After 30 years vegetation restoration, the area of forest lands 
significandy increased from 5% to 40% [24]. 

The Zhifanggou catchment has been an experimental site of the 
Institute of Soil and Water Conservation, Chinese Academy of 
Science (CAS) since 1973 [25]. The major agricultural land use 
type in LHR is slope cropland. Agricultural management in this 
region, including the major crop types grown, has not been 
changed significantly since the 1970s. The main crops grown in 
these sites were millet {Setaria italicd) and soybean [Glycine max) 
rotation, and no irrigation was provided in grown season (depend 
on rainfall). One crop was grown each year, and fertilizer was 
applied (mainly manure). After more than 30 years of compre- 
hensive management, the ecological environment of the catch- 
ment has been significantly improved [26]. Since late 1970s, slope 
cropland is replanted with shrubs and woods, mainly Robinia 
pseudoacacia L. (RP) and Caragana Korshimkii Kom (CK), to control soil 
erosion (see Table 1). Abandoned cropland was also generated 
during this period due to its extremely low productivity and long 





Vegetation types 


Age 


Canopy closure (%) 


Litter accumulation (t.ha ^) 


Undergrowth Vegetation* 


Species diversity indices 


Robinia pseudoacacia L 


30 


58 


20.5 


Lespedeza dahurica - Stipa bungeana 


6.5 


Caragana Korshinsl<ii Kom 


30 


50 


13.3 


Achillea capillaries, Stipa bungeana 


3.9 



*means the main vegetation in forest/shrub land. 
cloi:10.1371/journal.pone.0099657.t001 
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SOC content (g.kg-1) -p^ ^^ntent (g.kg-1) 



1 2 3 4 5 6 7 8 .2 .4 .6 .8 




POC content (g.kg-1) LOG content (g.kg-l) 
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Figure 2. Distribution of soil organic carbon (SOC, A), total nitrogen (TN, B), particulate organic carbon (POC, C), and labile organic 
carbon (LOC, D) contents of different land used types in soil depth of 0-200 cm. The error bars are the standard errors. 
doi:1 0.1 371 /journal.pone.0099657.g002 



distance from farmers' residences [27-28]. Despite wild grasslands 
and shrub lands were usually found on steep slopes, these sites 
were used for firewood collection as well. So the wild vegetation 
was of limited coverage or even barren for long periods. In 1999, 
most slope lands were closed for vegetation restoration under the 
GTGP [29]. 

Soil sampling 

In September 2012, based on land use history, 30 year old 
Robinia psendoacacia (RP), Caragana Korshinskii Kom (CK), abandoned 
land (AB) and slope cropland (SC) in the Zhifanggou catchment 
were selected. Three 30 mx20 m plots were established for each 
land use type. All sites were located on the same physiographical 
units with same slope aspects, same elevation of 1250 m and a 
spatial distance of 1200 m. 

Soil samples were taken at several soil depths using a soil auger 
(diameter 5 cm) from 10 points within "S" shape at each plot (0- 
10 cm, 10-20 cm, 20-30 cm, 30-40 cm, 40-50 cm, 50-60 cm, 
60-70 cm, 70-80 cm, 80-90 cm, 90-100 cm, 100-120 cm, 120- 
140 cm, 140-160 cm, 160-180 cm, and 180-200 cm). Then after 
removing the litter layer, ten soil samples at each depth of each 
plot were mixed to make one sample. Samples were collected at 
least 80 cm away from the trees. AH samples were sieved through a 
2 mm screen, and roots and other debris were removed. Soil 
samples were air-dried and stored at room temperature for the 



determination of soil chemical properties. A ring tube was used to 
determine the bulk density in each soil depth. 

Laboratory analysis 

SOC content (g.kg~') and TN content (g.kg~') were deter- 
mined usingK2Cr207 oxidation method and Kjeldhal method, 
respectively [30]. 

To determine POC content, 25 g soil was dispersed with 
100 mL of 5 g L~' sodium hexametaphosphate before being. 
Then, the mixed soil solution was shaken for 1 h at high speed on 
an end-to-end shaker and screened by a 0.053 mm sieve with 
several deionized water rinses. The soil remained on the sieve was 
backwashed into a pre-weighed aluminum box and dried at 60°C 
for 24 h, then it was grounded for analysis of C [31]. 

Soil labile organic carbon (LOC) was measured following the 
method described in Graeme et al. [32]. A 2-6 g air dried soil 
sample was put into a 50 mL centrifuge tube, and 25 mL of 
333 mmolL KMn04 solution was added before being shaken 
with a rate of 120 rpm for 1 h, and centrifuged for 5 min with a 
rate of 5,000 xg. The upper clear solution was transferred, and 
diluted by 250 times, and then the absorbance at 565 nm 
wavelength was determined. The absorbances at values 565 nm 
with different KMn04 concentrations were also determined for 
preparation of standard curve, which was used for the determi- 
nation of the KMn04 concentrations. Difference between the 



PLOS ONE I www.plosone.org 



3 



June 2014 | Volume 9 | Issue 6 | e99657 



SOC Pools in Deep Soil 





c 

O - 


^ 


4_ 












3 - 


a> 




o 














2- 






i 








o 


1 - 


o 




CO 






0- 




0.8 -] 








0.6- 














-•— • 








0 


0.4- 




















0.2- 


O 




O 




CL 






0.0- 



IRP 
ICK 

]AB 



a 

-Jil 




0-10 1040 40-100 100-200 

Soil depth (cm) 



0) 

o 
c 

0) 

I 

73 




0.2 



o 

O 




0-10 10-40 40-100 100-200 

Soil depth (cm) 



Figure 3. Differences in soil organic carbon (SOC, A), total nitrogen (TN, B), particulate organic carbon (POC, C), labile organic 
carbon (LOC, D) contents between SC and RP, CK or AB (RP/CK/AB - SC). Error bars are the standard errors. Different lowercase letters 
indicate significant difference among different land use types within same soil layer (P<0.05). The same for Fig 4 
doi:l 0.1 371/journal.pone.0099657.g003 



amounts of KMn04 added and remained was used to calculate 
labile C concentration in the soil sample. 

Calculation of SOC (TN) stocks, SR of SOC (TN, POC, and 
LOC) and CMI 

SOC density (SOCD, TND) represents the total SOC (TN) 
storage of overall certain sampling depth. SOCD(TND) of 
different samphng depths were calculated: 

SOCD(TND) = Csoc.TN xpxHx(l -5/100) x 10"' (1) 

where SOCD (TND) is the density (Mg-ha"') of SOC (TN) and 
Csoc.TN is the content (g'kg~') of SOC (TN). p is the bulk density 
(g-cm~^), H is the soil horizon thickness (cm), and 5 is the fraction 
(%) of gravels >2 mm in size in soil. Because the soil gravel size of 
loess in China is mosdy below 2 mm, this fraction was assumed to 
be 0 [33]. 

SR values (0-10 cm: 10-40 cm, 0-10 cm: 40-100 cm and 0- 
10 cm: 100-200 cm) were calculated from the contents of SOC, 
TN, POC and LOC following the method in Franzluebbers 
(2002). 

CMI values were calculated using following procedures: 
Firsdy, a C pool index (CPI) was calculated: 



CPI- 



sample total organic C (g/kg) 
reference sample total C (g/kg) 



(2) 



LI = 



L in each sampled soil 
L in the reference soil 



(3) 



where reference soil is SC soU, and L was calculated from the C 
lability: 



L = 



content of labile C 
content of non — labile C 



(4) 



At last, CMI was calculated: 



CMI = CPI X LI X 100 



(5) 



Statistical analyses 

All statistical analyses were carried out with SPSS 17.0. Analysis 
of variance (ANOVA) and Duncan's Multiple Range Test 
(DMRT) at 5% level of significance were used to compare the 
difference in contents and/or stocks of SOC, TN, POC, LOC, 
CMI, and SR among different land use types or soil depths. A 
sample linear-regression analysis was used to estimate the 
relationships between carbon stocks with CMI or SR values. 



where reference sample is SC soil. Then, a lability index (LI) was 
calculated: 
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SOC Stocks (Mg.ha-^) TN stocks (Mg.ha^) 




POC stocks (Mg.ha"^) Lqc stocks (Mg.ha" ^) 

0.0 0.5 1.0 1.5 2.0 2.5 3.0 o.O 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 




Figure 4. Stocks of soil organic carbon (SOC, A), total nitrogen (TN, B), particulate organic carbon (POC, C), labile organic carbon 
(LOC, D) of different land use types. The error bars are the standard errors. 
doi:1 0.1 371 /journal.pone.0099657.g004 



Results 

Changes in contents of SOC, TN, POC and LOC 

The contents of SOC, TN, POC and LOC responded 
differently as the change of soil depth (Fig. 2). In all land use 
types, contents of SOC, TN, POC and LOC in top soil (0-10 cm) 
were 3.26-7.86 g.kg"', 0.39-0.72 g.kg"', 0.65-1.31 g.kg"' and 
0.76-1.07 g.kg ', respectively, which were significantly higher 
than other soil layers (P<0.05). The contents of SOC, TN, POC 
and LOC decreased significantly in soil depth of 10-40 cm while 
the decreases trended to be flatter in subsoil (40-100 cm). 
Additionally, the differences in contents of SOC, TN, POC and 
LOC in deep subsoil (100-200 cm) were negligible (P<0.05). 

The differences in contents of SOC, TN, POC and LOC 
between three forest/shrub types (RP, CK and AB) and SC are 
shown in Fig. 3. The differences in SOC, TN, POC and LOC of 
RP and SC in soil depths of 0-10 cm and 100—200 cm were 
significantiy higher than that between other land use types and SC 
(P<0.05). The differences in SOC and TN of RP were 33.78% 
and 45.97% larger than that of CK, and 54.13% and 67.28% 
larger than that of AB in soil depth of 0-10 cm (P<0.05), while the 
differences in POC and LOC were 32.8%, 54.0% higher than that 
of CK, and 23.3% and 45.0% higher than that of AB (P<0.05). 



Moreover, the differences in SOC, TN, POC and LOC of RP 
were 25.05-85.29% higher than that of CK, and 61.78-90.70% 
higher than that of AB in soil depth of 100-200 cm. Additionally, 
significant differences in SOC, TN, POC and LOC contents were 
observed between RP and CK in soil depths of 10-40 cm and 40- 
100 cm, but there was no difference between CK and AB (P< 
0.05). 

Changes and distribution of SOC, TN, POC and LOC 
stocks 

SOC, TN, POC and LOC stocks of RP, CK and AB were 
higher than SC in aU soil profiles (Fig 4). The SOC, TN, POC and 
LOC stocks of RP were significantly increased (P<0.05), which 
were 0.43-5.8 Mg.ha"', 0.25-4.70 Mg.ha"', 0.44-9.14 Mg.ha"' 
and 1.49-11.38 Mg.ha"' higher than that of SC in soil layers of 
0-10, 10-40, 40-100 and 100-200 cm, respectively. Moreover, 
the stocks of SOC, TN, POC and LOC in soil layer of 100- 
200 cm of RP were higher than that of CK and AB by 15.4- 
32.1% and 21.8-43.1%, respectively (P<0.0,5). 

The SOC, TN, POC and LOC stocks responded differendy as 
the change of soil depth (Fig. 5). Although the distribution of SOC, 
TN, POC and LOC stocks in soil depths of 0-10 cm and 10- 
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Figure 5. Distribution ratios of soil organic carbon (SOC, A), total nitrogen (TN, B), particulate organic carbon (POC, C), labile 
organic carbon (LOC, D) in soil depth of 0-200 cm under different land use types. 

doi:1 0.1 371 /journal.pone.0099657.g005 



40 cm accounted for the majority, 26.36-34.06% and 21.08- 
36.62% were distributed in soil layers of 40—100 cm and 100- 
200 cm, respectively. Among four land use types, the highest 
proportion of SOC, TN, POC and LOC stocks were found in RP, 
while the lowest were in soil depths of SC in 0-10 cm and 100- 
200 cm of SC. The proportion of SOC, TN, POC and LOC 
stocks under RP were higher than SC by 4.68%, 7.32%, 4.65% 
and 5.96% respectively in soil depth of 0-10 cm soil depth, 
whereas by 5.90%, 9.78%, 6.30% and 10.06% was higher in soil 
depth of 100-200 cm soil depth respectively. 
Change in SR and CMI values 

Responses of SR in different land use types to change of soil 
depth were different (Fig 6). The SR values of SOC, TN and LOC 
differed significantiy among different soil depths (P<0.05), while 
the SR values of LOC diflfered only between 0-10:10-40 cm, 0- 
10:40-100 cm and 0-10: 100-200 cm. Among four land use types, 
the SR values of SOC, TN, POC and LOC of RP were the 



highest, but that of SC were the lowest in each soil depth(P<0.05). 
The SR values of SOC, TN, POC and LOC were in a decreasing 
order of CK>AB>SC. The SR values difllered significantiy 
between CK or AB with SC (P<0.05), while there was no 
significant difference between CK and AB. Additionally, the ratios 
of SR values of SOC, TN, POC and LOC in the surface layer (0- 
10 cm) to that in layer of 10—40 cm were >2.0 in most cases. 

The CMI values were significantly affected by land use types. In 
our study, the CMI values were in a decreasing order of RP> 
CK>AB>SC in four soil profiles and CMI values were 
significantiy enhanced by RP compared with SC (Fig 7). Averaged 
CMI values of RP, CK, and AB were 40.60%, 50.54%, 37.81%, 
and 14.1% higher than that of SC in soil layers of 0-10 cm, 10— 
40 cm, 40-100 cm and 100-200 cm. 

Regression equations to assess CMI/SR values of TN, POC, 
and LOC (Y) were showed in Table 2. There was a significant 
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0-10:10-40 0-10: 40-100 0-10: 100-200 

Soil depth (cm) 



Figure 6. Comparison of stratification ratio of soil organic carbon 
C), labile organic carbon (LOC, D) under different land use types. 

soil depths within same land use type while the different lowercase letters 
soil depth. The error bars are the standard errors. 
doi:10.1371/journal.pone.0099657.g006 

positive correlation between CMI/SR values of TN, POC, and 
LOC with SOC stocks in surface soil and deep soU. 

Discussion 

SOC, TN, POC and LOC contents and SOC and TN stocks 

Vegetation can greatly influence soil quality, C and N cycling, 
and regional socioeconomic development [34—35]. It is also 
reported that converting cropland into land with perennial 
vegetation would increase the SOC content [36]. Our results 
showed that land use type and soil depth significantly affected the 
contents of SOC and TN (Fig. 2). The conclusion that both land 
use type and soil depth are important factors influencing the soil 
carbon and nitrogen distribution was consistent with previous 
studies [35,37]. We also observed that the lowest SOC, TN, POC 
and LOC contents were found in slope cropland (Fig 4), which 
essentially agree with a previous study [38], indicating that the 
conversion of slope cropland to vegetation improves the C and N 
contents. A possible reason is that the lower residue input into the 
soil in slope cropland leads to lower SOC and TN contents. 




0-10:10-40 0-10: 40-100 0-10: 100-200 

Soil depth (cm) 



(SOC, A), total nitrogen (TN, B), particulate organic carbon (POC, 

Different uppercase letters indicate significant difference among different 
indicate significant difference among different land use types within same 



Additionally, our results showed diat SOC, TN, POC and LOC 
contents of RP were greater than that of CK and AB (Fig 4). It 
infers that the effects of RP on soil C and N play a significant role 
in land use and ecosystem management. The conclusion was 
consistent with Qiu et al [39], who reported that RP has potential 
to improve SOC content in the loessial guUy region of the Loess 
Plateau and the improvements are greater in long-term than 
middle-term. 

Recently it was reported that the depth of sampling is an 
important factor for the measurement of change in SOC stocks 
[40], and land use could influence subsoil C pools [41]. We found 
that SOC, TN, POC and LOC stocks of RP, CK, and AB were 
higher than SC for different soil profiles, especially in depths of 
40-100 cm and 100-200 cm (Fig. 4 and 5). It is demonstrated that 
converting slope cropland into woodland and shrubland not only 
affects SOC and TN stocks in surface soil, but also largely 
influences that in deep soil. The result was consistent with Wang et 
al [42], who reported that deep layer (50-200 cm) SOC stocks 
were equivalent to approximately 25% of that in the shallow layer 
(0-50) in HiUy Loess Plateau. This is mainly due to the fact that 
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0-10 10-40 40-100 100-200 

Soil depth (cm) 



Figure 7. Carbon management index (CMI) values of different 
land use types at different soil depths. The error bars are the 
standard errors. Different lowercase letters indicate significant differ- 
ence among different land use types within same soil depth. 
doi:1 0.1 371 /journal.pone.0099657.g007 

SOC input into subsoil is largely afTected by plant roots and root 
exudates, dissolved organic matter and bioturbation. In addition, 
most important factors leading to protection of SOC in subsoil 
include the spatial separation of SOM, microorganisms and 
extracellular enzyme activity related to the heterogeneity of C 
input [43]. As a result, stabihzed SOC in subsoil is horizontally 
stratified. 

Stratification ratios of SOC, TN, POC, and LOC 

According to Franzluebbers [7], SOC SR values >2 in 
degraded conditions is uncommon, and the SR values of SOC 



are generally low and seldom reach 2.0. And SR values of soil 
organic C and N pools with value of >2 would be an indicator 
that soil quahty might be improved [7] . In our study, the most of 
SR in SOC, TN, POC and LOC was more than 2 after convert 
slope cropland to forest or shrub land (Fig. 6). This means soil 
quality was improved in these afforested soils without disturbance. 
Greater C stratification ratios could be related to the fact that, 
during soil recovery by re-vegetation or land abandonment, soil 
was undisturbed thus reducing oxidation and favoring soil C [44] . 
The result was consistent with Sa et al [45]. Similar results were 
also reported by Moreno et al [46] and Franzluebbers [7], who 
reported that stratification of SOC occurs over time when soil 
tillage and disturbance is stopped and it is usually greater in 
undisturbed soils than in disturbed soils. In addition, the 
stratification may increase with time, and SOC, TN, POC and 
LOC contents are still aggrading but have not reached soil C 
saturation yet. That is the reason why SR values of SOC, TN, 
POC and LOC under CK and AB were higher but no significant 
difierences were observed compared with SC (Fig. 6). Sa et al [45] 
concluded that the SOC pool stabilization may be attained in 
about 40 years after long-term no-tillage adoption. 

Carbon management index 

CMI value was calculated to obtain indications of the C 
dynamics of the system and provide an integrated measure for 
quantity and quality of SOC [15]. Soils with higher CMI values 
are considered as better managed [47]. We found that CMI values 
were significantly enhanced by RP forest compared with CK, AB 
and SC in both surface soil and subsoil and deep soil (Fig. 7). Soil 
management under RP plot was more appropriate to improve the 
SOC status than other land use types. Similar result was reported 
by Qiu et al [39], who illustrated that RP forest has significantiy 
increased SOC, total nitrogen, ratio of carbon to nitrogen and 
ratio of carbon to phosphorus compared to other vegetation types. 
Our result showed that there were significant positive correlations 
between SOC stocks and CMI/SR in both surface soil and deep 
soil (Table 2). These findings showed that SR values of SOC, TN, 



Table 2. Regression equations among SOC stocks and CMI/SR for different soil layers. 





Axis 




Soil depth (cm) 


Equations 




Significant level 


X, CMI/SR* 


Y, SOC stocks" 










CMI 


SOC 


0-10 


Y = 53.5-I-23.20X 


0.91 


P = 0.048 






10-40 


Y=- 146-1-25.81 X 


0.97 


P = 0.027 






40-100 


Y = 24.2-l-8.60X 


0.93 


P = 0.023 






100-200 


Y = 77.5-I-4.05X 


0.93 


P = 0.046 


TN 


SOC 


0-10:10-40 


Y=1.12-I-0.08X 


0.93 


P = 0.047 






0-10:40-100 


Y=1.26-I-0.07X 


0.96 


P = 0.032 






0-10:100-200 


Y=1.67-I-0.07X 


0.98 


P = 0.017 


POC 


SOC 


0-10:10-40 


Y=1.18-I-0.06X 


0.97 


P = 0.023 






0-10:40-100 


Y=1.55-I-0.04X 


0.93 


P = 0.041 






0-10:100-200 


Y=1.74+0.04X 


0.92 


P = 0.047 


LOC 


SOC 


0-10:10-40 


Y = 1 .30-I-0.08X 


0.96 


P = 0.031 






0-10:40-100 


Y=1.95-I-0.05X 


0.97 


P = 0.025 






0-10:100-200 


Y = 2.53-I-0.05X 


0.98 


P = 0.013 



^CMI = carbon management index, SR = stratification ration, TN = SR of total nitrogen, POC = SR of particulate organic carbon, LOC = SR of labile organic carbon. 
''For the Y-axis, the SOC stocks (0-10 cm, 40-100 cm, and 100-200 cm) were used to analyze correlations between SOC stocks and SR values of TN, POC, and LOC (0- 
10:10-40, 0-10:40-100, 0-10:100-200). 
doi:l 0.1 371/journal.pone.0099657.t002 
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POM, and LOG, and CMI are suitable indicators for evaluating 
soil quality and C changes induced by GTGP in surface soil and 
deep soil. 

Conclusion 

In this study, the SOC, TN, FOG and LOG contents of RP, 
CK and AB in soil layer of 100-200 cm were higher than SG, 
especially for RP plot. Although die SOG, TN, POG and LOG 
stocks in soil layer of 100-200 cm were lower, there was more than 
27.38-36.62%, 25.10-32.91%, 2L59-31.69% and 21.08-26.83% 
of SOG, TN, POC and LOG stocks were distributed in 100- 
200 cm soil depth under RP, CK and AB. Meanwhile, the SR of 
SOG, TN, POG and LOG in the surface to lower depth ratio (i.e., 
0—10:10-40 cm) was >2.0 in most of case. And SR and as well 
GMI values were significandy enhanced by RP compared with SG 
in deep soil (100-200 cm) (P<0.05). Indicating that soil quality 
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